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Recent experimental results in direct dark matter detection may be interpreted in terms of a dark 
matter particle of mass around 10 GeV/c^. We show that the required scenario can be realized with 
a new dark matter particle charged under an extra abelian gauge boson Z' that couples to quarks 
but not leptons. This is possible provided the Z' gauge boson is very light, around 10-20 GeV/c^ 
in mass, and the gauge coupling constant is small, a' ~ 10~^. Such scenarios are not constrained 
by accelerator data. 



I. INTRODUCTION 

Many astrophysical observations support the existence 
of cold dark matter (CDM) in our universe, but its prop- 
erties, such as its mass, spin, and interactions, are still 
largely unknown. Several decades of experimental work 
have been spent trying to detect dark matter particles 
coming to Earth (direct detection), or more generally 
products of dark matter reactions in outer space (indirect 
detection). In these dark matter searches, one must dif- 
ferentiate signals due to dark matter from signals of other 
origin. In direct searches, a paramount telltale sign is a 
regular variation of the particle counting rate with a pe- 
riod of one sidereal year: an annual modulation m . Since 
1998, the DAMA experiment has been detecting such a 
modulation [2], with a current significance of more than 
8(7 [3]. Other direct detection experiments have placed 
strong constraints on the interpretation of the DAMA 
modulation in terms of dark matter particles (see e.g. 
for a recent analysis). 

A possibility that is not clearly excluded is a weakly- 
interacting massive particle of mass around a few 
GeV/c^, as proposed by Gondolo and Gelmini [S]. This 
is a natural mass range to explain the cosmic density of 
cold dark matter by means of thermal production in the 
early-universe plasma [5]. Recently, this possibility has 
become quite enticing thanks to reports of excess counts 
of unknown origin in the CoGcNT [7 and CRESST \8, 
experiments, and the hint that the CoGeNT events are 
modulated in a way similar to DAMA's Although 
these results are in tension with those of CDMS [10], 
XENONIO [n], and XENONIOO [12], cahbration uncer- 
tainties in the energy scale allow for a common interpre- 
tation of all the current measurements in terms of a ~7 
GeV/c^ dark matter particle with a ~10~'*'^ cm^ scatter- 
ing cross section with nucleons jl3j . 
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Particle physics models that include a light dark mat- 
ter particle with the properties just discussed must be 
extensions of the Standard Model (SM). The minimal su- 
persymmetric extension (MSSM), very light neutralinos 
might be possible at the edge of the available parameter 
space (see [13], but see also [TS] for opposite opinions). 
Other extensions of the SM may offer less tight possi- 
bilities (see e.g. [H]). In these and other models |17) . 
the interaction of the dark matter particle with quarks is 
(dominantly) mediated by the exchange of a light Higgs 
boson. Models with extra abelian gauge bosons Z' offer 
other possibilities for GeV dark matter. In these models, 
the coupling between the dark matter and ordinary mat- 
ter could be achieved through kinetic mixing of the Z' 
with the SM photon and Z boson p[S] , through exchange 
of extra fermions ^19J , or through the exchange of the Z' 
boson itself [ini HI]- 

In this paper, we examine Z' models in which the cou- 
pling of dark and ordinary matter is achieved through the 
exchange of the Z' boson itself. It turns out that in these 
models the DAMA/CoGeNT region is compatible with a 
correct dark matter density when the dark matter anni- 
hilates in the early universe through a Z' resonance. This 
requires ([55J[55) a relic density calculation more care- 
ful than previous ones. We find that the appropriate Z' 
boson should be light, ~ 10-20 GeV/c^, compatible with 
the supersymmetric U(1)b model in j21| but in contrast 
to the work of [50] who require a ~150 GeV/c^ Z' boson 
to describe the CDF dijet anomaly. 

II. EXTRA U(1)' gauge BOSON AS A 
MEDIATOR 

We consider an extra U(l)' gauge group under which 
the CDM particle and some standard model particles are 
charged. The U(l)' symmetry is spontaneously broken 
according to a nonzero vacuum expectation value (vev) 
of a Higgs field (j>. Then the massive U(l)' gauge boson 
becomes a mediator between the CDM sector and the 
SM sector. 

The couplings of the Z' boson to SM particles must 
be chosen wisely. If the SM leptons are charged under 
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U(l)' in a universal way, there is a strong constraint from 
LEP-II and Tevatron [2] on the ratio mz' / g' between 
the mass m^/ of the extra gauge boson and the couphng 
constant g' of the extra U(l) gauge group. For example, 
for U(1)=U(1)b_l, the LEP-II bound is approximately 
mz'/g' ^ 6 TeV/c^. The constraint can be relaxed if 
electrons are neutral under U(l)'. If only muons, taus 
and the CDM particles are charged under U(l)', one can 
construct a viable leptophillic dark matter model that 
accounts for some excess events observed in indirect dark 
matter searches |25j , but the required dark matter mass is 
much heavier than the few GeV/c^ we seek here. If only 
quarks and the CDM are charged (leptophobic Z'), the 
LEP-II/Tevatron constraints are even more relaxed [26l - 
128] ■ and as we show in this paper, it is possible to have 
a Z' model with a light dark matter particle. Finally, 
the Z' boson can couple to quarks and leptons indirectly 
through kinetic mixing with the SM photon and Z boson. 
Although such kinetic mixing is strongly constrained by 
experiments [571 [55], it can lead to viable models for dark 
matter interactions [TH [29] . 

Here we present two Z' models with a viable light 
CDM candidate: in one, the CDM particle is a scalar 
boson, in the other, it is a Dirac fermion. We realize 
the coupling to ordinary matter needed for direct CDM 
detection through the exchange of a leptophobic Z' (cou- 
plings of the Z' with the r lepton are in principle possible 
too). Although we impose no Z'-^ or Z'-Z kinetic mix- 
ing, a small amount of kinetic mixing could be generated 
radiatively. However, it turns out that in our scenario 
the gauge coupling constant g' is small (~10~^) and any 
radiatively-generated kinetic mixing is highly suppressed. 



A. Scalar CDM with Z' mediator 

Here we consider a scalar CDM particle X charged 
under the new U(l)' gauge group. We add the following 
lagrangian to the standard model, 

4,,!,, = D,X^ D^X ~ m\X^X - ^{X^X f 



X^XH^H - ^q^^q^X^X 



'^<f>UH^H-\z'^^z"^r 



(1) 



We also add a U(l)' term to all standard-model covariant 



derivatives Df^, 



Dl^-^Q'g'Z' 



(2) 



Here g' is the U(l)' gauge coupling constant and Q' is the 
U(l)' charge of the field on which acts. We will also 
use a' = g'^ /{Att). We assume that the U(l)' gauge bo- 
son couples to the SM fermions vectorially, and that the 



gauge anomaly is canceled by new fields charged under 
U(l)' [211 [30]. 

In the above lagrangian, H is the SM Higgs boson, 
and (p is the U(l)' Higgs boson whose vev breaks U(l)'. 
In order to realize (0) 7^ and (H) ^ 0, at least 
and must be negative. The Xhc/, term stabilizes the 
vacuum, and the negative squared masses must be tuned 
in correspondence with the size of the Xh4> term. 

This type of a model has been studied by Wise et 
al. [3D]. In their model, U(l)' is a gauged baryon sym- 
metry U(l)^, and Q'x is fixed by the Yukawa coupling 
that allows non-SM charged particles to decay. We will 
later discuss the case as a concrete example. 

For the stability of the CDM particle X on cosniolog- 
ical timescales, the field X should not acquire a nonzero 
vev {X), because a nonzero {X) induces the trilinear cou- 
pling \}ix{X)XH'^ H , which allows for decay processes 
such as X — > H^H. Further terms that cause the X 
particles to decay, such as Xcfr'^^/'^*' ^ arise from one- 
loop and non-renormalizable corrections. One can im- 
pose that the stability of the X particle at the renormal- 
izable level be means of the conditions 



Q'x ^ ±2Q' , iQ' 



(3) 



but one cannot completely forbid X to decay through 
higher-order couplings. We assume that m\ > is sat- 
isfied and that higher-order unsafe couplings are small 
enough to guarantee a life time for the X particles com- 
parable or greater than the age of the universe. In our 
numerical calculations, we take Q'x = 1, Q'q = 1/3 for 
quarks, and = for leptons, which is typical of a 
U(1)b coupled to baryon number. 

The annihilation cross section for XX* ff through 
Z' exchange is 



(7 



XX^ff 



S^Q'^Q'fa'^l3l3'{2E^ +m)) 
(4S2-m|,)2- 



J'J^ Z' 



(4) 



Here Q'j is the U(l)' charge of SM particle, /, with mass 
TO/, E is the X OT f center-of-mass energy, and (3 = 

— rn^x/ E'^ and /3' = y^l — m^^/E'^ are the center-of- 
mass velocities of each X and / in the center of mass 
frame. In the low-velocity limit, the leading term is a 
p-wave. The decay width of the Z' boson is given by 



mz 




12m z 



(5) 



In principle, the CDM particle X can annihilate to quarks 
not only through Z' exchange, but also through the Higgs 
bosons H and 0. For simplicity we neglect the Higgs 
boson contributions under the assumption that either the 
Higgs bosons are heavy or their couplings Ax0 and Xhx 
are small. 
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In addition, for mx > rnz', the annihilation into Z' 
pairs is allowed, with cross section 



'4 „/2 



'XX- 



16 



■Z'Z' " 

-8z^ 



w 



32 - 24z2 + 5^4 



4z2 + z4 ^ 4^2 

1 + (w + 



In 



i + {v-wy 



(6) 



Here z = mz'/mx, v = p/mx, and w — k/mx, where 
E and p are the center-of-mass energy and momentum 
of the initial X particles and k = yjE"^ — m?^, is the 
momentum of the final Z' particles. 



B. Dirac fermion CDM with Z' mediator 

Here we discuss the case of light CDM being a Dirac 
fermion charged under U(l)'. The SM is augmented by 
the following lagrangian terms. 



fermion 



V^x + a'Q'xf' - 'fnx) V'x 



T 



Here, ipx is a Dirac fermion with U(l)' charge Q'^. A 
global symmetry ipx — > e^^i^x can be enforced after the 
U(l)' symmetry breaking, so that ^x is guaranteed to be 
stable. However, higher-order non-renormalizable terms 
could generally break the global symmetry. Such terms 
are for example udd'ip^/ A? if U(l)' = U(l)^ and Q'^ — 1, 
and N^tpx/-^'^ if a right-handed neutrino Nn is added 
and charged under U(l)' with Q'^ = 1 and Q'x„ — —1/3. 
We assume that the cut-off scale A is large enough that 
the X particles are cosmologically stable. 

The annihilation cross section for XX ff through 
Z' exchange is 



'XX^ff 



E^ 



/3 (4i;2 



(8) 



In the non-relativistic limit, the s-wave contribution dom- 
inates. As for scalar CDM, we assume Higgs exchange 
contributions are suppressed either because the Higgs 
masses are large or because the X-Higgs couplings are 
small. 

The Z' width Tz' is given by Eq. (Sj) with the m|, — 
4m^ in the X contribution replaced By to^, + 2m^ . 

As for the scalar CDM case, for mx > mz' there is 
an extra contribution from annihilations into pairs of Z' 
bosons with cross section 



'XX^Z'Z' 



E^ V 
6-2z^ + z^ + 12v^ + 



(2 + z2)2 



2vw{l 



(2 + 2w2 - ^2)2 - 4wu; 

(9) 



4w\ l + (w + w)2"l 
In 



III. DIRECT DETECTION 

Since we assume that the U(l)' gauge boson couples to 
the SM fcrmions vectorially, the direct detection rate has 
a spin- independent component. Since the interaction is 
described by vector-current operators in the microscopic 
theory, the effective coupling of the Z' with protons, neu- 
trons and nuclei can be obtained using the conservation 
of the U(l)' charge. For a nucleus N of mass number 
A and electric charge Z, one has the effective lagrangian 
term 



Q'^Z'^N^'^N, 



where 



(10) 



(11) 



Q'x^ZQ'^ + {A-Z)Q'^ 
is the U(l)' charge of the nucleus, and 

Q'r,^Q'u + 2Q'^, Q'p = '2Q'u + Q'd (12) 

(7) are the U(l)' charges of the neutron and proton, respec- 
tively. For our choice of Q'^ — 1/3, we have Q'^ — A. 



f 

The non-relativistic limit of the spin-independent cross 
section for direct detection then follows as 



167ra'2 ,2„,2 / mxniN 

crxN = 4 — QxWn , 

m%, ynix + mis! 



(13) 



where m^v is the mass of the nucleus. 

Eq. (13) directly constrains the value of mz'/g' once 
mx and axp are determined in direct dark matter detec- 
tion experiments. For example, for Q'^ = 1 and ^ 1, 
the DAMA/CoGeNT region around mx ^ 1 GeV and 
c^xp ^ 10""*" cm2 leads to mz'/g' - 1 TeV. 

For simplicity, and as benchmark for our discussion, we 
take the DAMA/CoGeNT region outlined in [T3], shown 
in orange on the mx-crxp plane in Fig. IT] Fig. 1(a) 
corresponds to a Dirac fermion X, Fig. 1(b) to a scalar 
X (the direct detection constraints are identical in the 
two panels). Other analyses of the CoGeNT exponen- 
tial excess (e.g. [U[71[3T]) recover different regions in the 
iTix-o'xp plane, mostly to the right of the orange region 
shown, some compatible and some incompatible with the 
DAM A/LIBRA modulation region. 

Fig. [T] also shows the best current bounds from nega- 
tive dark matter searches (the excluded region in blue): 
CRESST at lower masses (from [4]) and XENONIO at 
larger masses. Two curves are shown for the XENONIO 
bound, one from [T^], the other from [i]. They reflect dif- 
ferent assumptions on the light detection efficiency near 
the threshold of the detector: the assumption with the 
higher detector sensitivity excludes the DAMA/CoGeNT 
region, the other does not. 
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IV. RELIC DENSITY 

The thermal density of the CDM particles X is given 
by the Boltzmann equation, 



Y decay width (a) CDM fermion 



dn 



3Hn = -{(T^nnV){n'^ 



"eg J 



(14) 



where n is the X number density and rieq is its value in 
thermal equilibrium. 

To compute the relic density, we use the procedure 
in [35' as implemented in DarkSUSY For this pur- 
pose, we introduced into DarkSUSY the invariant annihi- 
lation rate W = SE'pCTann, where CTann = J^f'^XXf^ff + 

cxxi^z'Z' is the total XX* or XX annihilation cross 
section given above. 

We impose that the computed cosmic density of X par- 
ticles Qxh'^ (in units of 1.8783 x 10"^^ kg/m^) equals the 
observed value of the cold dark matter density ilch^ = 
0.1123±0.0035 The thermal relic density depends on 
a', mz' and mx- If we fit the DAMA/CoGeNT region, 
the resulting parameters a' and mz' lead to a thermal 
density that is too large unless the annihilation is close, 
but not too close, to the resonance at mx ~ mz'/2. 

Contour lines of Vlx = in the mx-cxp plane are 
shown in Fig. [Tjfor several values of rriz' (the error bars 
on flch^ are within the thickness of the lines drawn). 
The parameter a' changes along each line. Below each 
line, one has fix > ^c- The thick red and purple con- 
tours correspond to mz' =12 GeV/c^ and 20 GeV/c^, 
respectively. Each contour shows a dip at mx = mz' /2 
due to the annihilation through the Z' resonance. As a 
function of mx, the resonance dip is highly asymmetric, 
being wider at mx < mz' /2. This is the correct behav- 
ior expected from the finite-temperature momentum dis- 
tribution of particles X during annihilation in the early 
universe (32]. 

We see that the fix — f^c contour lines sweep the 
DAMA/CoGeNT region for Z' masses in the range ^ 1 
to - 20 GeV/c^, touching the DAMA/COGeNT region 
on the left at the lowest mz' and on the right at high- 
est mz'- Fig. [2] gives a better visualization of the range 
of masses mz' and coupling constants a' that fit the 
DAMA/CoGeNT region. 

Notice that a heavy Z' with mz' ^150 GeV/c^, such 
as in suggested explanations of the CDF Wjj anomaly, 
has trouble matching the DAMA/CoGeNT region. If 
such a heavy Z' couples universally to quarks, fix in 
the DAMA/CoGeNT region would be too high, as seen 
by the location of the 150-GeV/c^ dashed line in Fig. [l] 
A correct X density may be obtained with non-universal 
couplings to quarks ^ Q'ujQ'd' summarily assumed 
in [,13], but such non- universal couplings may be very 
difficult to implement in a viable model without violating 
constraints from, for example, fiavor changing neutral 
currents. 




10 12 



mx [GeV/ci 




mx [GeV/ci 

FIG. 1: Contour lines of Qx = flc for several values of the 
Z' boson mass niz' ■ On each contour, the cosmic density of 
particles X (fermions in panel (a), and scalars in panel (b)) 
equals the cosmic density of cold dark matter. Also shown 
are the DAMA/CoGeNT region (in orange), direct detection 
constraints (in blue), and accelerator constraints (in yellow). 



ACCELERATOR BOUNDS 



As discussed above, direct detection requires mz' j g' ~ 
1 TeV/c^, while the fix — flc constraint leads to 
mx --10-20 GeV/c^ It follows that g' ~ lO'^, which 
is small but not unreasonably small. This region of 
small g' and small mz' is hard to reach in accelera- 
tor experiments, specifically because by assuming a lep- 
tophobic Z' we have avoided otherwise strong experi- 
mental constraints from LEP-II and the Tevatron. For 
U(1)'=U(1)b, the strongest bounds come from the invisi- 
ble and hadronic decay widths of the T meson 26, 28 , 34]. 
The region excluded by these bounds is shown in yellow 
in Fig. [T] the edges of the yellow region corresponding 
to the Z' masses plotted in the figure, namely mz' — 12 
GeV/c^ (lower edge) and mz' = 20 GeV/c^ (upper edge). 
The constraint from the invisible T width }34j is stronger 
at small mx , that from the hadronic width [2S] at larger 
mx- Clearly, accelerator bounds have no effect on our 
scenarios. 
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FIG. 2: U(l)' gauge boson mass m^/ and coupling constant 
a' that can explain the DAMA/CoGeNT region with a light 
cold dark matter scalar (in purple) or Dirac fermion (in red). 
Regions on the right have mx on the "left" of the resonance 
{mx < mz' /2); regions on the left have mx > mz' /2. 



imental results, i.e. have mass mx ~ 7 GeV/c^, scatter- 
ing cross section with nucleons axp ~ 10"^*^ cm^, and 
cosmic density equal to the cosmic density of cold dark 
matter. In one model the dark matter particle is a scalar, 
in the other a Dirac fermion. Both models assume that 
the interaction of the dark matter particles vifith ordinary 
matter occur through the exchange of a new leptophobic 
gauge boson Z' , and that the dark matter is produced 
thermally in the early universe. 

We find viable scenarios in which the Z' boson is light, 
with mass mz' ^ 10-20 GeV/c^, and gauge couphng con- 
stant g' ^ 0.02, or a' ~ 10~^, which is smaller than the 
Standard Model gauge coupling constants but not un- 
reasonably small. The small values of mz' and g' make 
accelerator constraints on our leptophobic Z' models in- 
effective. 
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